Current induced switching of magnetic domains to a perpendicular configuration 
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In a ferromagnet-normal-metal-ferromagnet trilayer, a current flowing perpendicularly to the layers 
creates a torque on the magnetic moments of the ferromagnets. When one of the contacts is super- 
conducting, the torque not only favors parallel or antiparallel alignment of the magnetic moments, 
as is the case for two normal contacts, but can also favor a configuration where the two moments are 
perpendicular. In addition, whereas the conductance for parallel and antiparallel magnetic moments 
is the same, signalling the absence of giant magnetoresistance in the usual sense, the conductance 
is greater in the perpendicular configuration. Thus, a negative magnetoconductance is predicted, 
in contrast with the usual giant magnetoresistance. 

PACS numbers: 75.70.Pa., 75.30.ds, 73.40.-c, 75.70.-i, 74.80Fp 
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A system made of stacks of alternating ferromagnetic 
and non-magnetic metal layers shows Giant Magneto- 
Resistance (GMR) When two consecutive magnetic 
layers have their magnetic moments aligned, the conduc- 
tance is much bigger than when they are anti-aligned. 
A simple way to think about this effect is in term of 
two separate (i.e., incoherent) currents for the majority 
and minority electrons (electrons with spin parallel or an- 
tiparallel to the magnetic moment, respectively), and to 
view the ferromagnetic layers as spin filters that have dif- 
ferent conductances for majority and minority electrons 
g. In the configuration where the magnetic moments 
are aligned, the majority electrons are well transmitted 
by both magnetic layers while the minority electrons are 
(mostly) reflected. When the moments are anti-aligned, 
the majority spin direction of one layer is the minority 
spin direction of its neighbor so that all the electrons are 
reflected. The role of the magnetic field is to align the 
initially anti-aligned magnetic moments, thus giving rise 
to an increase of the conductance. 

When the magnetic moments of the layers make an 
angle 9 different from or tt, the simple "two-current" 
picture does not hold anymore; instead a description in 
terms of a coherent superposition of spin up and spin 
down is needed. While the component of the spin flux 
in the direction of the magnetic moment m of a ferro- 
magnetic layer is conserved when an electric current is 
passed through the layer (since the fluxes of majority 
and minority spins are conserved individually), the spin 
flux perpendicular to m does not have to be conserved, 
as it depends on the coherence between majority and mi- 
nority electrons. As first pointed out by Slonczewski || 
and Berger Q , the consequence of such a change of the 
magnetic moment carried by the current is that the cur- 
rent exerts a torque on the moments of the ferromagnets. 
This so-called "spin-transfer" torque vanishes at the two 
angles 9 = and 9 = n, where the two current model ap- 
plies. Experimentally, one studies a ferromagnet-normal- 
metal-ferromagnet (FNF) trilayer where one of the mag- 
netic moments is held fixed, e.g. by using a thick fer- 



romagnetic layer |3]-||, while the other one is free to ro- 
tate. The current direction determines which of the two 
angles 9 = and 9 = tt is stable; switching the current 
direction reverses the relative orientation of the magnets, 
and thus changes the conductance. This signature of the 
current-induced spin-transfer torque, which was observed 
experimentally in Refs |^,^], provides a mechanism for a 
current controlled magnetic memory element. The spin- 
transfer torque can also be used for other applications 
Q , including the excitation of spin- waves [|l0| . 

If one of the two external contacts of an FNF trilayer is 
replaced by a superconductor (S), see Fig. |], the picture 
changes drastically. For voltages smaller than the super- 
conducting gap, transport through the system occurs by 
Andreev reflection: an electron impinging on the S in- 
terface is reflected as a hole with opposite spin, adding 
a Cooper pair to the superconducting condensate. Each 
layer is therefore traversed twice, once by a majority elec- 
tron (or hole) and once by a minority hole (or electron). 
Hence, the 9 = and 9 = tt configurations have the same 
conductance, the difference between them only being the 
order in which the spin filtering occurs JTlj ] . 

What happens for angles 9 other than or tt? This 
question has a remarkable answer. While the presence of 
Andreev reflection suppresses the "usual" GMR at 9 = 
and 7t (see previous paragraph), we have found that it 
leads to a richer variety of effects for other 9, where quan- 
tum coherence between spin up and spin down is cru- 
cial. With the S contact, not only parallel or antiparallel 
alignment of the magnetic moments plays a special role, 
but also the perpendicular configuration: (1) As with 
normal contacts, passing a current through the FNFS 
system exerts a torque on the moments of the F layers; 
however, the possible stable configurations are not only 
9 = and tt, but also 9 = tt/2 and 3tt/2. (2) The conduc- 
tance g is a 7r-periodic function of 9 with a maximum for 
9 = tt/2. These two results were derived for the exper- 
imentally relevant case when the ferromagnetic moment 
rfib neighboring the superconductor is held fixed, while 
the other moment rh a is allowed to vary, see Fig. |l|, and 
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FIG. 1. Schematic drawing of the system. Normal metal 
layers are Ni, N2, and N3, the ferromagnetic layers with mo- 
ments rh a and mt in the y-z plane are denoted with F a and 
-F&, and the superconductor is denoted with S. The layers N 
are metallic spacers, and fictitious ideal leads 1, 2, 3, 4 and 5 
are added for technical convenience. 



assuming that the spin-relaxation length is larger than 
the system size, which is a reasonable assumption for a 
thin trilayer system. By switching the current direction, 
one can switch the orientation of the free magnetic mo- 
ment rh a from parallel to perpendicular to rhb- Since 
g(0) ^ g(7r/2), such a switch can be observed through a 
change in the conductance, similar to what was observed 
111 Refs. |,| for the case of normal contacts. A magnetic 
field will align the moments, bringing the system from 
9 ~ n/2 to 9 — (instead of 9 = it to 9 = 0), so that one 
should observe a positive magnetoresistance, in contrast 
to the standard GMR. 

The FNF trilayer under consideration is shown in Fig. 
[]]. The current flows in the ^-direction, perpendicularly 
to the layers. For technical reasons, we have added thin 
ideal (non-magnetic) spacers labeled I, . . . , 5 between the 
F and N layers. Denoting the spin current in layer i by 
Ji (i = 1, . . . , 5), the torques fb and f a on the magnetic 
moments rh a and fhb are given by 
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(1) 



Equation (Q) expresses that, unlike the probability cur- 
rent, the spin current J is not conserved by the ferro- 
magnetic layers. (The component of J in the direction of 
the magnetization is conserved, though.) The difference 
in spin currents is transfered to the F layers, as a torque 
acting on their magnetic moments. 

We first give an intuitive explanation of how the pres- 
ence of the superconductor allows for the additional sta- 
ble configurations at 9 = tt/2 and 9 = 3tt/2, and then 
present the results of a more rigorous calculation. For the 
intuitive explanation we make the simplifying assump- 
tion that spins scattered from an F layer with moment 
fh have their magnetic moment pointing parallel (or an- 
tiparallel) to fh. (This assumption is valid if majority 
and minority electrons are transmitted incoherently; It 
is not necessary for the stability of 9 = 7r/2 or 37r/2, as 
shown in our calculations below.) With this assumption, 
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FIG. 2. Left panel: signs of the torques (small arrows) 
for different orientations of the moments fh a and fht, (long 
arrows). Upper (lower) panel: right (left) moving electrons. 
The torques lie in the plane spanned by fh a and fh . Right 
panel: schematic of the geometrical argument leading to 
f a — in the perpendicular configuration. 

the explanation proceeds as follows (see right panel of 
Fig ||). (i) For voltages below the superconducting gap, 
only Cooper pairs can enter S. Hence, no spin current can 
flow into S and one must have J4 = J5 = 0. Therefore, 
by Eq. (|]), fb = J3. (ii) Since the spin flux in the direc- 
tion of rhb is conserved, there can be no component of the 
torque fb directed along rhb- Hence J \ = J3 is perpen- 
dicular to rhb- (iii) Finally, for an unpolarized electron 
current entering the trilayer, any spin current in layer 1 
must be carried by electrons scattered from F a , so that 
Ji is parallel to rh a . The torque f a is then given by the 
component of Ji perpendicular to rh a . Hence, 



\Tb\ \ COS I 



(2) 



f a thus vanishes at 9 = n/2 and 37r/2 yielding the new 
possibility of stable configurations at these angles. 

To make the above picture more precise, we use the 
scattering approach |l4|,|l5|]. The trilayer is bounded in 
the y and z directions, so that its transverse degrees of 
freedom are quantized, giving iV c h propagating modes 
("channels") at the Fermi level. One has N c ^ ~ A/Xp, A 
being the cross section of the system and Xp the Fermi 
wave length. Expanding the electronic wave function in 
these modes, we can describe the system in terms of the 
2A r c [ 1 -component vectors v[/^' l ) i (- R ) ; which is the projec- 
tion of the wave function onto the left (right) going modes 
in region i for the the electrons (holes) (i — 1, ... ,5). The 
2iV c h components of 'J account for the spin and channel 
degrees of freedom. In the scattering approach, each layer 
is characterized by 2iV c h x 2N c b reflection matrices r% and 
r[ and transmission matrices t[, U (the label i = 1 for F a , 
i = 3 for Fb, and i — 2, 4 for the two normal spacers), 
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At the NS interface electrons are reflected as holes, 
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To find the electrical and spin currents I and Ji for each 
layer i, we need to calculate the generalized 27V c h x 2A r c h 



2 



scattering matrices S^}^^ that give the amplitudes of 
left (right) moving electrons (holes) in region i in terms 
of the amplitudes of the incoming electron (hole) in the 
normal electrode (labeled with layer index i = 1), 



e{h)L(R) 



ge(h)e(h)^e(h)L 
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(5) 



MR) 



can be expressed in 



Using Eq.(||), the matrices S. 

terms of r^, r[, t[ and fj. (This was done in Ref. [|l3j for 
an FNF trilayer with two normal-metal contacts.) 

Following the derivation of the Landauer formula for 
the conductance one can now calculate the electri- 
cal current / fl7|| and the spin current Ji for an applied 
voltage V across the system as 

dl _ 2e 2 
dV ~ 9 ~ ~ 

dJi 
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where a = (a x , a y , a z ) is the vector of Pauli matrices. 

It is important to notice that the scattering matrices 
need to be calculated for each individual layer only sep- 
arately; solution of Eqs. (§) - (||) then describes how 
the information from the individual layers is combined 
to give properties of the multilayer system. For the fer- 
romagnetic layers, the precise form of the reflection and 
transmission matrices has to be determined from a mi- 
croscopic model, see e.g. Ref. |16 . For the normal metal 



spacers, we make use of the polar decomposition MA 
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= Is 
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where i — 2,4, I2 is the 2x2 unit matrix in the spin 
grading, U{ u^, Uj, and v[ are A^h x iV c h unitary matrices, 
and Ti is a diagonal matrix containing the eigenvalues 
of tit] on the diagonal. We consider an ensemble of tri- 
layers for which the metallic spacers are good diffusive 
metals (with different impurity configurations for differ- 
ent members of the ensemble), or for which the interfaces 
between the different material are rough, so that the dis- 
order mixes the different transverse modes in an isotropic 
way. In that case one can describe the ensemble by taking 
the unitary matrices u, u' , v, and v' uniformly distributed 
in the unitary group |14|j . We calculate an average over 
the matrices u, u' , v, and v' to find the ensemble average 
of dl/dV and dJi/dV. Since, fluctuations are of rela- 
tive size l/iV c h, the average is sufficient to characterize a 
single sample when iV c h 1. JTJ] (In the experiments, 
typically N ch ~ 10 3 .) 

After the ensemble average, the results only depend on 
four parameters for each ferromagnetic layer, [ fi"3"[ 

T iT = N^tttiAv T ^ = N ^nAv 
Tin = N- h hr t n tl , R tU = N- h Hr rqrj x , (8) 
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FIG. 3. Conductance (top) and torque per unit current 
(bottom) as a function of 9 for a Co-Cu-Co trilayer system. 
dr/dl is expressed in units of h/Aive 2 , g in units of N c he 2 /h. 
We took T aT = T fet = 0.68 and T al = Tm = 0.29. These 
numerical values were obtained from Ref. |16J. The conduc- 
tance of both normal spacers is taken equal to 2A r c i 1 e 2 /h, cor- 
responding to a situation where the main effect of disorder 
in the spacers is to randomize directions, not to backscat- 
ter electrons. The transmission/reflection coefficients Tfj and 
Rfl have been set to zero. The solid (dashed) curve shows 
r a (r(,) in the presence of the S contact. For comparison, the 
dot-dashed curve shows Tb = r a for two normal contacts jL3|. 



on the conductances of the normal spacers, and on the 
angle 0. In Eq. (0), the arrows refer to the majority and 
minority spin directions. The coefficients T-jq and i?iq de- 
scribe the coherence of transmission and/or reflection of 
minority and majority spins. (The quantity N c ^ — N^R-^ 
is known as the "mixing conductance" JOl.) The as- 
sumption of "incoherent" transmission and reflection of 
majority and minority that we made in the intuitive ar- 
gument above, amounts to setting T-jq and i?-q to zero. 
This is a fair assumption for T-|q, since majority and mi- 
nority electrons pick up different and uncorrelated phase 
shifts upon transmission through the ferromagnets, but 
it does not have to be the case for when the reflec- 
tion at the interface is instantaneous. For "perfect" spin 
filters, where all minority spins are reflected and all ma- 
jority spins are transmitted, T-jq and i?-|q are both zero. 

The resulting expressions for the conductance and the 
torques are rather lengthy, and will not be reported here. 
Nevertheless, they allow us to draw general conclusions 
about the direction and magnitude of the torques for ar- 
bitrary values of the system parameters. First, in agree- 
ment with the simple argument given above, we find that 
the torques f a and t% lie in the plane spanned by the mo- 
ments rh a and rhb [i.e., parallel to (m a x rhb) x rh a and 
(fhb x rh a ) x rhb respectively]. Second, the magnitudes of 
the torques f a and 77, are related as 
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irrespective of the conductances of the normal layers. 
This equation replaces Eq.(||) in case of non vanishing 
Rati and T -|-j_. In Fig. |[ we have indicated the direc- 
tions of T a and Tf, for various relative orientations of fn a 
and rfib and for the two possible directions of the current. 

In Fig. ||, the torques f a and t% and the conductance 
g are shown versus 8, for realistic choices of the scat- 
tering parameters of the Co/Cu/Co trilayer used in the 
experiment of Ref. |||| , as they can be obtained through 
ab-initio calculations, see Rcf. p"6| . The torque and con- 
ductance for normal contacts with the same scattering 
parameters are also shown for comparison. Note that the 
presence of the superconducting contact does not affect 
the order of magnitude of the torques, or the sensitiv- 
ity of g to 9. Therefore the critical current necessary to 
switch the layers' magnetic moments should be in the 
same range as for N contacts (lent ~ 10 9 A/cm~ 2 in the 
point contact geometry of Ref. || and ~ 10 7 A/cnU 2 for 
the pillar geometry of Ref. H). This would allow, at 
least for the pillar geometry, the use of bias voltages be- 
low the superconductor gap. The main difference from 
the case of normal metal contacts is the period of the 
^-dependence of f Q and g, which changes from 2tt to it. 
The ^-dependence of the torque j% on the layer adjacent 
to S is similar to what is found for the case of two normal 
leads. Although Fig. || is for a special choice of the pa- 
rameters, we have verified that these conclusions do not 
depend on the detailed choice of scattering parameters. 

In the limit when the conductance g^ of the normal 
spacer connecting F a and Fj, is small, and when the mixed 
transmission and reflection probabilities Tj-j. and R^i are 
<C 1, simple expressions can be obtained for the torque 
and the conductance, 
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Note that, in this limit, f a , t%, and g do not depend on 
the detailed properties of layer b but only on the direction 
of its magnetic moment rfib. 

The predicted absence of a GMR in the usual sense 
[g(9) is equal for 8 — and tt] requires some clarification, 
as GMR has been observed in multilayers with supercon- 
ducting contacts fli| , in contrast to our predictions and 
those of Ref. (T^| . One possible reason for the observed 
GMR is spin relaxation between the ferromagnetic layers 
and between the ferromagnetic layer and the supercon- 
ductor. Especially for multilayer systems, spin relaxation 
can not be neglected. Another cause can be that ferro- 
magnets not only serve as spin filters, but they also cause 
scattering from minority electrons into majority electrons 
and vice versa Hj. Finally, for ballistic FNFS layers, 
quantum interference effects that are not considered here 
can be important, thus providing another mechanism for 



the conventional GMR In each of these cases, the 

conductance of the multilayer will be different for 8 = 
and 8 = tt, although the underlying reasons for this GMR 
are more subtle than in the case of normal metal contacts. 

To conclude, we would like to emphasize that the sta- 
bility of the perpendicular configuration is an intriguing 
feature; for example, it allows for the construction of a 
nanomagnet resonator with a current and magnetic-field 
controlled frequency: when a current is passed through 
the system such that the torque f a keeps the moment 
m a in the plane perpendicular to rhb, an applied mag- 
netic field parallel to fhb will make rh a precess in its plane 
with a frequency controled by the magnetic field. 
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